PERIODIC SOLUTIONS IN AN ARRAY OF COUPLED 
FITZHUGH-NAGUMO CELLS 
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o 

» I , Abstract. We analyze the dynamics of an array of A''^ identical cells coupled 

C^ ' in the shape of a torus. Each cell is a 2-diniensional ordinary differential equa- 

tion of FitzHugh-Nagumo type and the total system is Zjy x Zjv— symmetric. 
The possible patterns of oscillation, compatible with the symmetry, are de- 
Qs , scribed. The types of patterns that effectively arise through Hopf bifurcation 

*vj ' are shown to depend on the signs of the coupling constants, under conditions 

' ensuring that the equations have only one equilibrium state. 

q 

'^ ■ 1. Introduction 

<^ . 

Hopf bifurcation has been intensively studied in equivariant dynamical systems 

in the recent years from both theoretical and applied points of view. Stability of 

equilibria, synchronization of periodic solution and in general oscillation patterns, 

- stability of the limit cycles that arise at the bifurcation point are among the phe- 

J^ I nomena whose analysis is related to the Hopf bifurcation in these systems. Periodic 

fv^ i solutions arising in systems with dihedral group symmetry were studied by Golu- 

l^ I bitsky et al. [10] and Swift [8], Bias and Rodrigues [2] dealt with the symmetric 

^D ' group, Sigrist [12] with the orthogonal group, to cite just a few of them. Bias et 

tH- I al. [T] studied periodic solutions in coupled cell systems with internal symmetries, 

f^ . while Bionne extended the analysis to Hopf bifurcation in equivariant dynamical 

CO I systems with wreath product ^3^ and direct product groups Jy . The general theory 

of patterns of oscillation arising in systems with abelian symmetry was developed 

by Filipsky and Golubitsky [5] . The dynamical behavior of 1-dimensional ordinary 

differential equations coupled in a square array, of arbitrary size (2iV)^, with the 

symmetry Z^v x Z^v, was studied by Gillis and Golubitsky [5]. 

C^ . In this paper we use a similar idea to that of [6] to describe arrays of N'^ cells 

where each cell is represented by a subsystem that is a 2-dimensional differential 

equation of FitzHugh-Nagumo type. We are interested in the periodic solutions 

arising at a first Hopf bifurcation from the fully synchronised equilibrium. To each 

equation in the array we add a coupling term that describes how each cell is affected 

by its neighbours. The coupling may be associative, when it tends to reduce the 

difference between consecutive cells, or dissociative, when differences are increased. 

For associative coupling we find, not surprisingly, bifurcation into a stable periodic 

solution where all the cells are synchronised with identical behaviour. 
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When the coupHng is dissociative in either one or both directions, the first Hopf 
bifurcation gives rise to rings of N fuhy synchronised cehs. All the rings oscillate 
with the same period, with a -^-period phase shift between rings. When there is 
one direction of associative coupling, the synchrony rings are organised along it. 
Dissociative coupling in both directions yields rings organised along the diagonal. 
The stability of these periodic solutions was studied numerically and were found to 
be unstable for small numbers of cells, stability starts to appear at iV > 11. 

For all types of coupling, there are further Hopf bifurcations, but these neces- 
sarily yield unstable solutions. 

This paper is organized as follows. The equations are presented in section [2] 
together with their symmetries. Details about the action of the symmetry group 
Zjv X Zjv are summarized in section |3l we identify the Z^v x Z^r— irreducible sub- 
spaces of R^^ ; the isotypic components; isotropy subgroups and their fixed point 
subspaces for this action. This allows, in section |4l the study of the Hopf bifurca- 
tion with symmetry Zat x Zn , applying the abelian Hopf bifurcation theorem [S] 
to identify the symmetries of the branch of small-amplitude peridic solutions that 
may bifurcate from equilibria. In section [S] we derive the explicit expression of the 
2N'^ eigenvectors and eigenvalues of the system linearized about the origin. Next, 
in section |6] we perform a detailed analysis on the Hopf bifurcation by setting a pa- 
rameter c to zero. In this case the FitzHugh-Nagumo equation reduce to a Van der 
Pol-like equation. Finally, in section [71 we characterize the bifurcation conditions 
for c > small. 

2. Dynamics of FitzHugh-Nagumo coupled in a torus and its 

symmetries 

The building-blocks of our square array are the following 2— dimensional ordinary 
differential equations of FitzHugh-Nagumo (FHN) type 

(--^N X ^x{a-x){x-l)-y = fi{x,y) 

ij^hx-cy = f2{x,y) 

where a,b,c ^ 0. Consider a system of A^^ such equations, coupled as a discrete 

torus: 

(2) 

Xa.,l3 = Xa,l3 (« - Xa,p) {Xa,li " 1) " J/q./S + 7(2^0, /3 " a;Q+l,^) -|- 5{Xa,p - Xa,p+l) 
ya,l3 = bXa^p - Cya,l3 

where 7 7^ (5 and l^a^A, l^/3^iV, with both a and /3 computed (mod N). 
When either 7 or J is negative, we say that the coupling is associative: the coupling 
term tends to reduce the difference to the neighbouring eel, otherwise we say the 
coupling is dissociative. We restrict ourselves to the case where A ^ 3 is prime. 

The first step in our analysis consists in describing the symmetries of ([2]). Our 
phase space is 

K^^' = {{xc.,fi,yo.^p) I 1 < a,/3 < TV, x„^^,2/a,^ e M} 

and (21) is equivariant under the cyclic permutation of the coUumns in the squared 
array: 

(3) 7i(a;i,,3, • ■ ■ , XN,i3; yi,p, ■ ■ ■ , yAT./j) = {x2,p, ■ ■ ■ , xn,i3, a:i,/3; ?/2,/3, ■ ■ ■ , yN,p, yi,p) 
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as well as under the cyclic permutation of the rows in the squared array: 

(4) 72(a;aa, ■ • • ,Xa,N;ya,l, ■ ■ ■ ,ya,N) = {Xa,2, ■ ■ ■ , Xa,N , Xa,l] ya,2, ■ ■ ■ .yaM.Vas)- 

Thus, the symmetry group of ^ is the group generated by 71 and 72 denoted 
Zjv X Zjv = (71, 72)- We will use the notation 7[ • 7I e Zat x Zjv as (r, s) = 7[ • 7|- 
Let's refer to the system ([2]) in an abbreviated form i — f{z), z — {xa,i3,ya,p) 
while A S K is a bifurcation parameter to be specified later. The compact Lie group 
Zjv X Zjv acts linearly on R^^ and / commutes with it (or is Z^ x Z^v—equi variant) . 

We start by recalling some definitions from [TU] adapted to our case. 

The isotropy subgroup Sz of Zjv x Zjv at a point z S R^' is defined to be 

T,z — {(r, s) G Zjv X Zn : (r, s) • z = z} . 
Moreover, the fixed point subspace of a subgroup S G Z^v x Z^v is 
Fix (E) = {z e R^^' : (r, s) • z = z, V(r, s) e e| 
and /(Fix(S)) CFix(S). 

Definition 1. Consider a group T acting linearly on R". Then 

(1) A subspace V C R" is said T— invariant, if a ■ v £ V , Vcr G F, Vw e V; 

(2) ^ subspace V C R" is said F— irreducible if it is T— invariant and if the 
only T~invariant subspaces of'V are {0} and'V. 

Definition 2. Suppose a group T acts on two vector spaces V and W. We say 

that V is T -isomorphic to W if there exists a linear isomorphism A : V — > W 
such that A{ax) = aA{x) for all a; S V. // V is not T -isomorphic to W we say 
that they are distinct representations ofT. 

3. The Zn x Zn action 

The action of F = Zn x Zn is identical in the Xa.p and the ya^p coordinates, 
i.e. F acts diagonally, ^{x,y) — (7a;, 7?/) in R^^ for x,y € R^ . Hence, instead of 
taking into account the whole set of {xa,i3,ya,^) G R^^ , we will partition it into 
two subspaces, R^' x {0} and {0} x R^^ namely x^^p G R^' and y^^p € M^^ 
The action of Zn x Zn on K^ has been studied by Gillis and Golubitsky in [6| for 
N = 2n, we adapt their results to our case. 

Let k = (fci, fc2) G ^^ and consider the subspace Vk C M.^ , where (xa^p) G Vk 
if and only if 

'2ni 



(5) Xa.p = Re zexp 



N ("'^) 



nN^ 



zeC, l^a,(3 ^N 



Proposition 1. Consider the action o/Zjv x Zjv on R^ given in ([3]) and (j4]) with 
N prime and let I be the set of indices k — (fci,fc2) listed in Table[l\ Then for 
k G / we have 

(1) dimVk = 2 except for k = 0, where dimVo = 1. 

(2) Each Vk defined in ([5]) is Zn x Zn— invariant and Zn x Zn— irreducible. 

(3) T/ie subspaces Vk are a// distinct Zn x Zjv representations. 

(4) T/ie group element (r, s) acis on Vk as a rotation: 

"27ri 



(r, s) • z = exp 



(r, s) • k 

N ^ ' ' 
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(5) The subspaces Vk verify ®kG/ ^k = R^ • 

(6) The non-trivial isotropy subgroups for "Ln x Zat on M^ are Zjv ('', s); ^^e 
subgroups generated by one element {r,s) 7^ (0,0). 

(7) //(r,s)^(0,0) then 



Fix I Zjv (?", s) ) = 2_^ "^k o,nd dimFix I Z^v (?', s) 



N. 



k-(r,s) = 

(mod N) 



Table 1. Types of Zat x Zjy— irreducible representations in Vk G 
R^ , where Vk is the subspace of R^ corresponding to k = 

(fci,fc2) G/CZ^. 



Type 


dim (Vk) 


k 


Restrictions 


(1) 


1 


(0,0) 




(2) 


2 


(0,fc2) 


li^k2^{N-l)/2 


(3) 


2 


(fci,0) 


l^ki^{N^l)/2 


(4) 


2 


(fcl,fcl) 


1 «; fci s; (A^ - 1) /2 


(5) 


2 


(fcl,fc2) 


l^fc2<fcl^A^-l 



Proof. The arguments given in [B] Lemma 3.1] with suitable adaptations show that 
statements [T]-[4l hold and also that the subspaces Vk with k ^ / are redundant. 
Since the Vk are all distinct irreducible representations, a calculation using[T] shows 
that X)ke/'^™'^k = N'^, establishing El 

For [6l note that since N is prime, the only non trivial subgroups of Zat x Zat 
are the cychc subgroups Zat (r, s) generated by (r, s) ^ (0, 0). Each (r, s) fixes the 
elements of V(j,j j,^) when (fci,A;2) = {N — s,r). Using |31 it follows that (r, s) ^ 
(0,0) fixes X — (xa.p) G V(fcj_fe2) with x ^ if and only if exp [^ (?", s) ■ (^1,^2)] = 
1 i.e. if and only if (r, s) ■ {ki, fc2) = (mod A^). Thus Fix (r, s) is the sum of all the 
subspaces Vj^^ ^.^j such that {r,s) ■ (fci,fc2) = (mod N), it remains to compute 
its dimension. 

Let {k[,k2) — {N — s,r), so that V/^,, j,, \ C Fix(r, s). Then for any {ki,k2) we 

^^ ^^ (mod N). Thus V(fe^,fe,) C Fix (r, s) if and only 



have (r, s) • (/ci, fc2) — det 



if det 



ki k2 
ki ^2 



ki ^2. 
(mod N), and {r,s) G ker 



/ci fc2 
"^1 ^^2 



(mod iV). This is equiv- 



alent to having (A:i,fc2) and {k'i,k'2) linearly dependent over Zat, i.e. {k'i,k'2) — 
m{ki,k2) (mod iV). Then Fix (r, s) — Emgz„V„(Ar_s,r)- Since half of the (fci, fc2) — 
m{N — s,r) have fc2 > ki, this expression adds two times the same subspace. 
Then, since for {ki, ^2) ^ (0, 0) we have dimV(;;j ^j^) = 2 we obtain dimFix (r, s) = 

l + i(7V-l)-2 = iV. ' D 

3.1. The Zat x Za^— isotypical components of M^^ . So far we have obtained 
the distinct Zat x Zat— invariant representations of M , by considering only the 
subspaces corresponding to the variable a;. If V C M^ is an irreducible subspace 
for the action of Zn x Zn then V x {0} and {0} x V are Zat x Zjv— isomorphic 
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irreducible subspaces of M^^ . We will use the notation V © V for the subspace 

Vx VcK^^'. 

Definition 3. Suppose a group T acts on R" and let V C M" be a T -irreducible 
subspace. The isotypic component of M" corresponding to V is the sum of all T- 
irreducible subspaces that are T -isomorphic to V. 

Once we have the Zjv x Z^r— ireducible representations on R^ , we can calculate 
the isotypic components of the representation of this group on M^^ . 

Table 2. Isotypic components Zk of the Zjv x Z^v— action on 
E^^ , where the Vk are defined in ([5]). 



dim(Zk) 


Zk 


Restrictions 


2 

4 
4 
4 
4 


V(o,o) © V(o,o) 
V(fc,,o)©V(fe,,o) 


1 ^ fci s^ (TV - f ) /2 
I ^ fca < /ci s^ iV - f 



Lemma 1. The isotypic components for the Zjv x Z^— representation in M^^ are 
of the form Zk = Vk ® Vk given in Table\^ 

Proof. Since T — Zjv x Z^v acts diagonally, 7(2;, y) = {jx,jy) in E^^ for x,y £ 
E^ , then the subspaces Vk®{0} and {O}0Vk are F— invariant and irreducible, by 
H of Proposition[TJ If k =^ k' then, by El we have that Vk' ® {0} and {0} ® Vk' are 
not F— isomorphic to either Vk ® {0} or {0} ® Vk- Therefore, the only isomorphic 
representations are Vk ® {0} and {0} ® Vk for the same k G / and the result 
follows. D 

4. Symmetries of generic oscillations patterns 

The main goal of this section is to characterize the symmetries of periodic solu- 
tions of ([2]), specially those that arise at Hopf bifurcations. 

Given a solution z{t) with period P of a Zjv x Z^v-equivariant differential equation 
i = /(z), a spatio-temporal symmetry of z(t) is a pair (cr, 0), with a G Z^r x 1,^ 
and 6* e E (mod P) ~ S^ such that a ■ z{t) = z{t + 9) for all t. The group of spatio- 
temporal symmetries of z{t) can be identified with a pair of subgroups H and K 
of Zat x Zat and a homomorphism : 7J — >■ S^ with kernel K , where H represents 
the spatial parts of the spatio-temporal symmetries of z{t), while K comprises the 
spatial symmetries of z{t), i.e. the symmetries that fix the solution pointwise. In 
order to get all the spatio-temporal symmetries for solutions of ^, we use the 
following result of Filipitsky and Golubitsky: 

Tiieorem 1 (abelian Hopf theorem [5]). In systems with abelian symmetry, gener- 
ically, Hopf bifurcation at a point Xq occurs with simple eigenvalues, and there 
exists a unique branch of small- amplitude periodic solutions emanating from Xq- 
Moreover, the spatio-temporal symmetries of the bifurcating periodic solutions are 
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H = TjXo O'^d K ~ keiv(^), where V is the centre subspace of the bifurcation at 
Xq and H acts H -simply on V . 

Recall that a subgroup H of Zat x Z^r acts H— simply on a subspace V if either 
V is the sum of two isomoprphic iJ-irreducible subspaces or V is iJ-irreducible but 
not absolutely irreducible. 

Proposition 2. Periodic solutions of ^ arising through Hopf bifurcation with 
simple eigenvalues at an equilibrium point Xq have the spatio-temporal symmetries 
of TableWi 



Table 3. Spatio-temporal (H) and spatial (K) symmetries of ^ 
that may arise through Hopf bifurcation at a point Xq, with two- 
dimensional centre subspace V. Here V^j- denotes the subspace 
Fix (Zn (k)) C R^' as in the proof of H of Proposition [H 



H 


set 


centre 


K 


restrictions 






containing Xq 


subspace V 




on k 




r 


Vo®Vo 


y C Vo ® Vo 


H 






r 


Vo®Vo 


F g Vk ® Vk 


I^N (k^) 


k=^0 




Z^(k) 


Vk^eVkAW 


1/ = Vo © Vo 


H 


k^^O 




Zw(k) 


Vk^eVuAJo} 


V<^Ve®Vi 


1 


k^^O 


ik^O 


ZAr(k) 


Vk^eVkAW 


VgV,®Ve 


H 


k^O 


£k = 


1 


K2^'\UVk®Vk 


^ C Vk ® Vk 


1 







Proof. The proof is a direct application of Theorem [l] using the information of 
Section [3l From assertions El and [71 of Proposition [l] the possibilities for H are 
1, Zjv {r,s) and Z^r x Zat. This yields the first coUumn in Table El The second 
collumn is obtained from the list of corresponding fixed-point subspaces. 

Let V be the centre subspace at Xq. Since the eigenvalues are simple, V is two- 
dimensional and is contained in one of the isotypic components. Then either H acts 
on V by nontrivial rotations and the action is irreducible but not absolutely irre- 
ducible, or H acts trivially on y C FixTf , hence V is the sum of two i?-irreducible 
components. In any of these cases H acts if-simply on V. The possibilities, listed 
in Lemma HI yield the third collumn of Table El The spatial symmetries are then 
obtained by checking whether V meets Fix(_ff). D 

Hopf bifurcation with simple eigenvalues is the generic situation for systems with 
abelian symmetry [5l Theorem 3.1]. In the next section we will show in Theorem |4l 
that this is indeed the case for Q and obtain explicit genericity conditions. 

A useful general tool for identifying periodic solutions whose existence is not 
guaranteed by the Equivariant Hopf Theorem, is the H mod K Theorem [9]. Al- 
though it has been shown in [9] that in general there may be periodic solutions 
with spatio-temporal symmetries predicted by the H mod K Theorem that cannot 
be obtained from Hopf bifurcation, this is not the case here. 
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5. Linear Stability 

In this section we study the stabihty of solutions of ([2]) lying in the full synchrony 
subspace V(o q) ffi V(o,o) C M^^ . For this we choose coordinates in R^^ by con- 
catenating the transposed coUumns of the matrix {xa^p,ya.i3), i-e. the coordinates 
are (Ci, . . . , Cn)'^ where 

(6) C'lS ^{xi^p,yi^p,X2,i3,y2,l3,---,XN,l3,yN,ls) ■ 

Let p e V(o,o) ffi ^(0,0) be a point with all coordinates {xa,i3,ya.j3) = {x*,y*)- In 
these coordinates, the linearization of ^ around p is given by the N x N block 
circulant matrix M given by 

'A B ... 
A B ... 
A B ... 



M 



B 







A 



where ^4 is an TV x A^ block circulant matrix and S is an iV x iV block diagonal 
matrix given by 



A 



D 


E 
















D 


E 
















D 


E 








B 



F 



















F 



















F 











£; ... D 

where the 2x2 matrices E and F are given by 



E = 



-7 




F = 





-5 0" 




OOF 



and D is obtained from the matrix of the derivative D(fi, /2) of ([T]) at (x*,?/*) as 
D = D{fi, /2) — E — F. In particular, if p is the origin we have 

\d -ll 



D = 



with 



d = —c 



-1 



Given a vector u e C'^ , we use the iV*'* roots of unity w'' = exp {2Trir/N) to 



define the vector Cl{r, v) e 



■^kN 



n{r,v) = [w,w''D,a;2'-w,...,w(^-i)'^t;] , s^r s$ iV- 1. 
The definition may be used recursively to define the vector S(r, s,v) — Cl (s, r2(r, v)) G 



as 



S(r, s, v) = \n{r, v), Lj'Q{r, v),uj^'Q{r, «),..., J'^-^^'n{r, v) 

Theorem 2. If Xr,s is an eigenvalue and t; G C" an eigenvector of D + uj^E + uj'^F, 
and if M is the linearization of ([2]) around p £ V(o,o) © ^(0,0) then Ar,s is an 
eigenvalue of AI with corresponding eigenvector 'E.{r, s,v). 

Proof. Let us first compute the eigenvalues of the matrix A. We have, for any 

weC2 

An{r,v) ^n{r,{D + uj''E)v) 
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or, in full, 



D E ... 
ODE ... 
D ^ ... 



E 



D 



{D + uj''E) V 
uj'' {D + uj^'E) V 
Lo^'- {D + uj-'E) V 



so, if {D + uj^E) V = XrV then An{r, v) — Xrfl{r, v). 

By applying the same algorithm we can calculate the eigenvalues and eigenvectors 
of the matrix M. Given u e C^^ compute 



Mn{s, u) = n (s, (A + uj'B) u) 



or, in full 



A B ... 
Q A B ... 
Q Q A B ... 



B 



A 



UJ u 



{A + Lu''B)u 
uj' {A + CO" B) u 
uj'^" {A + Lo'' B) u 

uj(^-^)^A + uj''B)'. 



To complete the proof we compute ME{r, s, v) = Mfl (s, 51(r, v)) as 

ME{r, s,v) ^n (s, {A + uj^B) n{r, v)) = Q (s, An{r, v)) + fl (s, uj^Bflir, v)) . 

Then, since i? is a block diagonal matrix, then Bil{r, v) — n{r, Fv) for any w e M^ 
and we get: 

ME{r, s,v) =n (s, n (r, (D + oj^'E) v)) + fl (s, fl{r, w^Fu)) 
= nls,n (r, (D + uj^'E + uj"F) v)) 

and thus 

M'~{r, s,v)=E (r, s, {D + lo'' E + uj"F) v) . 

It follows that if [D + lu'^E + lo^F) v — Xr,sV then A/S(r, s, v) = Ar.sS(r, s, v) as we 
had claimed. D 

5.1. Form of the eigenvalues. 

Theorem 3. // L is the linearisation of ([T]) around the origin and D = L — E — F 
then the eigenvalues of D + cu'^E + lo'^F are of the form 

1 



A(r,«)± - -[-(c + a)+7(l-w'') + <5(l-w^)] 



(7) 



± J\/[(^ ~ ") + ^(1 " ^'') + -^(l " '^')]^ " "^^ 



where \/~^ stands for the principal square root. Moreover, on the isotypic component 
'^(ki,k2)®'^{ki.k2}^ t^^ eigenvalues of M are Xi^kiM)^ "'^'^ their complex conjugates 

^{N~kuN-k2)±- 

Proof. It is straightforward to derive the explicit expression ([7|) of the eigenvalues of 
D + uj'^E + uj'^F; a direct calculation shows that, unless (r, s) = (0,0), the complex 
conjugate of X{r,s)+ is not A(r,s)-, but rather X(^N-r,N-s)+- 
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We claim that for any complex nmiiber ^ the real and imaginary parts oi^{ki , ^2 , C) 
lie in 'V(ki.k2)j from this it follows that the real and imaginary parts of the eigenvec- 



tors S(fci,fc2,w) lie in the isotypic component 'V(ki.k2 
k2), this will comp 
N-s)± correspond to \(kiM) 



''^(ktM)- Since V 



{kiM) - 

'y^ (M-ki,N-k2)j this '^ill complete the proof that the eigenvalues A/j.^ ,t2)± and 



\N-r,N-s)± correspond to Y^^iM) ® ^(kiM)- 



It remains to establish our claim. Using the expression ([5]) to write the coordinate 

Xa,i3 of S(fci,fc2,C) 'we obtain 



Xa,l3 



Cexp 
Cexp 



(a-l)fci + (/3-l)fc2 

'2TTi 



N 



(-fci - fc2) 



exp 



N 



(a,/3).k 



Its real and imaginary parts are of the form ([S]) for z = Cexp [2^ (— fci — fc2)] and 



as claimed. 



-zCexp[fi(-fci 



^2)], respectively, and therefore lie in 'V{ki,k2) ® ^ 



(fel,fe2)' 

D 



For the mode (r, s) = (0, 0) the expression ([7]) reduces to the eigenvalues X± of 
uncoupled equations ([T]), linearized about the origin, 



(8) 



A± 



-c-a± ^y{c + ay -46 



6. Bifurcation for c — 

In this section we look at the Hopf bifurcation in the case c = 0, regarding a as 
a bifurcation parameter. The bulk of the section consists of the proof of Theorem[5] 
below. Since in this case the only equilibrium is the origin, only the first two rows 
of Table [3] occur. 

Theorem 4. For generic j,6 and for c = 0, h ^ Q all the eigenvalues of the 
linearization of ([2]) around the origin have multiplicity 1. 

Proof. We can write the characteristic polynomial for L + (w'' — 1)E + {uj'^ — 1)F, 
where L is the linearization of ([T]) about the origin, as 

(9) /(A,r,s)= X^ + \[a--i{l-Lo'')-5{l-Lo')]+b 



l<r,s^N. 



We start by showing that if two of these polynomials have one root in common, 
then they are identical. 

Indeed, let 0i(A) and (t>2{X) be two polynomials of the form ^ and suppose they 
share one root, say X = p + iq, while the remaining roots are A = pi + iqi for 0i 
and A = P2 + *'?2 for 4'2- Since 5 7^ 0, then none of these roots is zero. Then we can 
write 

0j(A) = {X - {p + iqj) {X - {pj + iqj)) 

= X^-X{{p + p,) + i{q + q,)) + ip + iq) {p, + iq.j) 

and therefore {p + iq) (pi + iqi) = b — {p + iq) {p2 + 192), so, aa p + iq ^ 0, then 
{pi + iqi) = (P2 + i92) and therefore 01(A) — 02(A). Since this is vahd for any 
pair of polynomials of the family, it only remains to show that for generic 7, S the 
polynomials do not coincide. 

Two polynomials /(A, r, s) and /(A, f, s) of the form (j9]) coincide if and only if 



(10) 



7 (w'^ - uj'~') = 6 {uj' - uj') . 
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Thus, for (7, 6) outside a finite number of fines defined by (fTOj) aff the eigenvafues 
of the finearization of ([2]) around the origin have muitipiicity 1, as we wanted to 
show. D 

Theorem 5. For c = 0, b > and for any 7 and 5 with 7(5 ^ the origin is the 
only equilibrium of ([2]) . For each value of 7 and S there exists a* ^ such that 
for a ^ a^ the origin is asymptotically stable. The stability of the origin changes 
at a — a», where it undergoes a Hopf bifurcation with respect to the bifurcation 
parameter a, into a periodic solution. The spatial symmetries of the bifurcating 
solution and the values of a, are given in Table CI Moreover, if the coupling is 
associative, i.e, if both 7 < and 5 < 0, the bifurcating solution is stable and the 
bifurcation is subcritical. 

Table 4. Detaiis of Hopf bifurcation on the parameter a for The- 
orem [5] Solutions bifurcate at a = a* (where 6m — — ^ ) with 
spatial symmetry K and spatio-temporal symmetry Zjv x Zjv- 

sign(7) sign((5) a* K 

Ztv X Zjv 

+ - 7(l~C0s6lAr) ZAr(0,l) 

+ (5(1 -cos Sat) ^ ZAr(l,0) 

+ + (7 + (5)(l-COS0^) ZAr(i^,^) 



The first step is to determine the stability of the origin. To do this, we need 
estimates for the real part of the eigenvalues ([7]) . This is done in the next Lemma. 

Lemma 2. LetA{r,s) = -a + 7 (1 - u;'')-h(5 (1 - w"*). i^or c = 0, & > 0, 7^ 7^ and 
for all {r,s) we have Re A(r.s)- ^ ^Re ^(^js)- -^/Re A{r,s) ^ then Re X(r,s)+ ^ 
Re A{r,s), otherwise Re A(r.s)+ < 0. 

Proof of Lemma In order to evaluate the real and imaginary parts of eigenvalues 
•A(r,s)±: 'W6 need to rewrite equation ([7]) by getting rid of the square root. For this 
purpose, we use a well known result from elementary algebra; we have that if 
r] = ai + ibi, where ai and bi are real, &i ^ 0, then the real and imaginary parts of 



y/r] = a/oT^M&i are given by 



(11) ReV^=yM±^ ImV^^sgn(Mv/^. 

A direct application of ([7]) in Theorem [3] to the case c = yields, if ^^ ^ M and for 

£1 =±1 



1 / . /M2-46|+Re (^2-46) 

(12) ReA(,,,),, = - ReA + en'' ' ^ ' 



2 



nox T X 1 /t . ^ /|A2-46|-Re(A2-46) 

(13) Im A(,^,),^ ^ - I Im A + eiS\ ' ' 

with A = A{r, s) and S = sgn (im [A^]) 
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The statement for A(r.s)- follows immediately from (|12p. For A(r,s)+; note that, 
for 6 > and any 77 G C, we have 

\ri^ - 46| + Re {if - 4fe) s^ 2 (Re -qf 

with equality holding if and only if Re 77 = 0, when the expressions are identically 
zero. Hence, taking rj = A(r, s), we obtain from (J12p : 

Re A(^^s)+ s$ 2 ^^^ ^^"^^ *) + 1^^ ^'^''' *)') 
and the result follows. D 

The particular case of fully synchronised solutions in Theorem [S] is treated in 
the next Lemma. This case is simpler since the bifurcation takes place inside the 
subspace V(o,o) ®V(o,o)- 

Lemma 3. For c = 0, fe > and 7(5 7^ the origin is the only equilibrium of 
^ and at a = it undergoes a Hopf bifurcation, subcritical with respect to the 
bifurcation parameter a, to a fully synchronised periodic solution. If both 7 < and 
(5 < 0, the origin is asymptotically stable for a > and the bifurcating solution is 
stable. Otherwise, the periodic solution is unstable. 

Proof of Lemma O Inspection of ^ when c = shows that the only equilibrium 
is the origin. 

The restriction of ^ to the plane V(o^o) © V(o.o) obeys the uncoupled equations 
(m whose linearisation around the origin has eigenvalues given by ([8]). It follows 
that, within V/g^o) © ^(0,0): the origin is asymptotically stable for a > 0, unstable 
for a < 0. The linearisation has purely imaginary eigenvalues at a = 0. 

Consider the positive function (p{y) = exp(— 2y/6). Then, for the uncoupled 
equations ([T|) we get 

d d 

— {(p{y)x) + — {(p{y)y) = (-Sx^ + 2ax - a) ip{y) 

which is always negative if < a < 3. Hence, by Dulac's criterion, the system ([1]) 
cannot have any periodic solutions, and thus if there is a Hopf bifurcation at a = 
inside the plane V(o_o) ® V(o,o) it niust be subcritical. 

In order to show that indeed there is a Hopf bifurcation we apply the criterium 
of [HI Theorem 3.4.2] and evaluate 

-L05 — J XXX I Jxyy i 9xxy 'T' Qyyy 

H JP [fxy \Jxx + fyy) ~ 9xy \9xx + 9yy) ~ JxxQxx + Jyy9yy\i 

where f{x,y) = fi{x,Vby) - Vby, 9{x,y) = f2{\/bx,y) ~ Vbx and fxy denotes 



(0, 0), etc. Since, for a = c = 0, we have f{x, y) = —x + x and g{x, y) = 



d'f ^ 
dxdy^ 

this yields s* — — |. The Hopf bifurcation is not degenerate and the bifurcating 
periodic solution is stable within V(o,o) © V(o,o)- Since 9Re X(ofi)±/da = —1/2 the 
bifurcation is indeed subcritical with respect to the bifurcation parameter a. 

It remains to discuss the global stability, with respect to initial conditions outside 
V(o,o) © "^(0,0)- If 7 > then from the expression p2|) in the proof of Lemma [2] at 
a = we obtain Re A(i_o)+ > ^-nd the bifurcating periodic solution is unstable. 
A similar argument holds for 6 > 0. 
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If both 7 < and 6 <0, then, for a = 0, we get Re A{r, s) < for (r, s) ^ (0, 0). 
Hence, by Lemma [21 all the eigenvalues X(r.s)±, {r,s) ^ (0,0) have negative real 
parts and the bifurcating solution is stable. D 

Proof of Theorem O The case of associative coupling 7 < and 5 < having been 
treated in Lemma [Sj it remains to deal with the cases when either 7 or (5 is positive. 
Let r^, = ^^, £2 = ±1, c = 0, 6 > 0, On = ^^^^^, with sin6'Ar > 0, cos^at < 
and cos ^^v < cos ^^ for all a e Z. 

If 7 > 0, 5 < and a ^ 7 (1 — cos^at) — a,, then for aU (r, s) ^ ('"eai 0) we have 
Re A{r,s) ^ 0, with equality only if both a — a^, and {r,s) — {re^^O)- Using (fT2|) 
and (fT3|) we get 

A(re2,0) = i£27shi6'Ar; 



ImA(r,^ ,o)ei = 2 ( ^27 sin 9n + £1 \/7^ sin^ ^JV + 46 



with A(r_|_.o)+ — -^(r_.o)- and A(r_,o)+ = -^(r+.o)-- In addition ImA(r^.o)+ > 
ImA(r_,o)+ > for 6 > 0. 

The results of Golubitsky and Langford [7] are always applicable to the Hopf 
bifurcation for A(r , ,o)+j since there are no eigenvalues of the form fcA(r_|_,o)+ with 
k € N. For the smaller imaginary part there may be resonances when A(r, ,0)+ = 
kXtT-,0}- with k € N. Otherwise, if the other non-degeneracy conditions hold, 
there are two independent Hopf bifurcations at a = a* i.e. two separate solution 
branches that bifurcate at this point. The resonance condition may be rewritten as 

72sin6'Ar - ^''~^' b, ken, k^2. 
k 

The bifurcating solutions are stable if and only if the branches are subcritical. The 
eigenspace corresponding to these branches lies in V(^r-±,o) © ^(r±,o) <^ Fix(Z(0, 1)). 
In the case 7 > 0, J > 0, we have a* = (7 + 5) (1 — cosOm)- For a ^ a„ 
and for all (r, s) we have Re A(r, s) ^ 0, and hence Re \{r.s)± ^ with equality 
holding only when both a — a^ and A(r_,_,r±)+- The eigenspace in this case lies in 

V(,^,,^)©V(,^,,^)CFix(Zjv(^,^)). Then 

A{ri.^,r^^) = 1(^27 + £3(5) sin 6* AT 

A(r,,,r,3),, = 2 (^('^£2 , ^e3 ) + ei V^^ (^£2 : Te^ ) - 4^) , 

with Ei — ±1, i — {1, 2, 3}. Hence, ReA/^e ,re )ei — and 



2 ( (£27 + £3(5) sin6lAr + eiy (£27 + £3^)^ sin^ On + 46 



ImA(r,^,r,3)s 

with \ir r \- — Afr r , , • In addition ImA/^ „ v, > when £1 = +1 with 

l'-e2''-e3/ '.'e2i'E3)+ We2''e3/^1 

luY\r+,r+)+ > IniA(r^^,r,3)+ ^^"^ ('^2,£3) ^ (+!,+!) and 6 > 0. 

Hence there is a non-resonant Hopf bifurcation corresponding to \{rj^,rj^)+- As 
mentioned before, for the smaller imaginary parts there may be resonances when 
^{re ,1-5 )+ = fc'^(r+,r+)+ with A; G N. Otherwise, there are four independent Hopf 
bifurcations at a = a* if other non-degeneracy conditions hold; in this case four 
separate solution branches bifurcate at this point. The bifurcating solutions are 
stable if and only if the branches are subcritical. D 
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We have checked numerically the non-degeneracy condition for bifurcation of 
the non-resonant branch using the formulas of Golubitsky and Langford [7]. The 
criticality of the bifurcation branch seems to depend on N. For 7 > 0, ^ > 0, 
A^ ^ 11, the bifurcating solution branch seems to be always subcritical, and hence 
stable. For N — 3,5,7 it seems to be supercritical. If (5 ^ 0, and N ^ 11, the 
bifurcating branch seems to be subcritical for large values of 7 > 0, supercritical 
otherwise. 

7. Bifurcation for c > small 

In this section we extend the result of section |6] for bifiurcations at small positive 
values of c. We start with the case when both 7 and 6 are negative. 

Corollary 1. For small values of c, z/ 6 > 0, 7 < and S < the origin is an 
equilibrium of ([2]) and there is a neighbourhood of the origin containing no other 
equilibria. There exists a near such that for a ^ a the origin is asymptotically 
stable. The stability of the origin changes at a — a, where it undergoes a Hopf 
bifurcation, subcritical with respect to the bifurcation parameter a, into a stable 
periodic solution with spatial symmetries Zjv x Zjv. 

Proof. The eigenvalues of Df{0) are all non-zero at c = 0, as shown in the proof of 
Lemma [31 Hence Df{0) is non-singular and the implicit function theorem ensures 
that, for small values of c, there is a unique equilibrium close to the origin. From 
the symmetry it follows that this equilibrium is the origin. 

If both 7 < and 6 < Q, then it follows from the proof of Lemma [3] that for 
c = the purely imaginary eigenvalues at a = a* = are simple. Continuity of the 
eigenvalues ensures the persistence of the purely imaginary pair for c 7^ at nearby 
values of a . Hence the corresponding eigenvectors depend smoothly on c, and the 
non-degeneracy conditions persist for small values of c. D 

The cases when either 7 > or (5 > 0, are treated in the next proposition. 

Proposition 3. For small values of c > 0, if b > and j6 ^ the origin is an 
equilibrium of ^ and there is a neighbourhood of the origin containing no other 
equilibria. For a > a*, where a* has the value of Table^ the origin is asymptotically 
stable. For almost all values 0/7 7^ and 6 j^ 0, The stability of the origin changes 
at a = a < a,, , with a, near a*, where it undergoes a non-resonant Hopf bifurcation 
into a periodic solution having the spatial symmetries of Table[^ If the bifurcation is 
subcritical with respect to the bifurcation parameter a, then the bifurcating periodic 
solution is stable. 

Proof. For small values of c, the origin is locally the only equilibrium, by the argu- 
ments given in the proof of Corollary [1] In Lemma[4] below, we show that for small 
c > and a > a,:, all the eigenvalues of the linearisation have negative real parts. 
Hence, the origin is asymptotically stable. When a decreases from a* the real parts 
of some eigenvalues change their signs. It was shown in the proof of Theorem [5] 
that for c — 0, there are several pairs of purely imaginary eigenvalues at a = a*. 
In Lemmas [5] and [6] below, we show that, generically, for small c > 0, when a de- 
creases from the value a* of Table [H the first bifurcation at a = a < a* takes place 
when a single pair of eigenvalues crosses the imaginary axis at a non-resonant Hopf 
bifurcation. We also identify the pair of eigenvalues for which the first bifurcation 
takes place. D 
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Lemma 4. For small values of c > 0, if b > 0, ^S ^ 0, let a* have the value of 
Table\^ If a '^ a^, then, for all r,s, and fore i = ±1, we have IieX(^r,s)ei "^ ^' '^'^'^ 
the origin is asymptotically stable. 

Proof. For ei — ±1, the eigenvalues X(^r,s)ei have the form 



46. 



(14) 2A(,^,),, = A{r, s)-c + ei^{A{r, s) + cf 
Using (|lip and writing A{r, s) = x + iy, we have 

(15) (A + c)'^ = (c + x)'^ -y'^ + 2i{cy + xy). 
Then ReA(r.s)- ^ ReA(r ^j^, and ReA(,. ,,.)+ ^ if and only if 

Rey{A{r, s) + cf - 46 s$ (x + c) 

This never happens if c — x < 0, and in this case we also have 2ReA(r_s)_ ^ 0. If 
c — a; ^ let 

pi = (c- x)^ - 4a;c + y^ +46 and p2 = [(c + a;)^ - ?/ - 46] +4y^(c + a;)^ 

with p\-p2 = 16 [-x^c + (6 + 2c^)x^ - (2c6 + cy'^ + c^)x + c^b] . 
With this notation, ReA(r,s)+ ^ if and only if pi > and 

[{c + xf ^y"^ -Ahf + Ay^{c + xf ^ [{c - xf - Axe + y"^ + Ab]^ . 
We have the following cases: 

(1) if X < 0, c > then pi > and Pi — _P2 > and so ReA(r,s)+ < 0; 

2 

(2) if a; = 0, c > then pi ^ y"^ + Ah + c^ > Q and ^1—^1 ^ c^b > and 

16 
therefore ReA(r_s)+ < 0; 

(3) at a; = c we have p\—p2 = — 16c^y^ < 0, so {p\ — P2) changes sign for some 
3^^ 5 U \ X^ \ c. 

This completes the proof, since for a ^ a* we have x = ReA(r, s) ^ 0, as in the 
proof of Theorem [S] D 

Lemma 5. For small values of c > 0, i/ 6 > 0, 7 > 0, (5 < and if a ~ a* < 
is small, then all the eigenvalues of the linearisation of ^ around the origin have 
real parts smaller than the real part 0/ A(r, ,o)+ 

Proof. It was shown in the proof of Theorem [5] that at c = 0, a = a», the eigen- 
values Xi^c ,o)£iJ with £1 ~ ±lwith £2 = ili are purely imaginary, and all other 
eigenvalues have negative real parts. From the expression ((T4| it follows that 

2ReA(,^^ ,o).i = ReA(r,, , 0) - c + c^Rc^ {A{r,,,Q) + cf - Ab 
hence ReA(r^ .0)- < ReA(r^ ,o)+- Let a be the value of a for which the pair 
A(r+,o)+ — ^(j--.o)+ first crosses the imaginary axis. Since Rej4(re2,0) decreases 
with a, then a < a^,. The estimates above show that at a = a, the second pair 
A(r, ,0)- — A(r_^o)- still has negative real part. For small c, the other eigenvalues 
still have negative real parts at a, by continuity. D 

Lemma 6. For small values ofc>0, ifb>0, ^>0, S>0 and if a — a* < 

is small, then all the eigenvalues of the linearisation of ([2]) around the origin have 
real parts smaller than the real part 0/ A(r+.r+)+ 
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Proof. As in Lemma [5] we use ()14|) to show that ReA(r^ .r^ )- < ReA(r^ ^^e )+ &t 
c > 0. It remains to compare the real parts of the two pairs X(^r+,r+)+ = -^(r_,r_)+ 
and \{r+,r-)+ = A(r_,r+)+- To do this, we write A{r^^,rf;^) = x + iy, where a;,y € R 
and obtain conditions on x and y ensuring that the eigenvalue is purely imaginary. 
Then we evaluate these conditions on the expressions for x and y to obtain the 
result. 

From ([Mil and (fTS]) we get that ReA = if and only if 

{2{c-xf -Xf =X'^ + 4:{x + cfy'^ for X = {x + c)^ - y'^ - ib 
and this is equivalent to 

which may be rewritten as: 



^c-xf-{x + cf+4b = y^ ^''^''^^ 



(c - xy 



This may be solved for y to yield 



2 _ (b - cx){c - x)'^ _ ( b , 2 



(16) y'='- '^^ ^^ _^i (c-a;)^=V(a;) 

ex \cx J 

and note that, for c > and if c^ < 6, then ip{x) > for < a; < c, and in this 
interval ipix) is monotonically decreasing. 

Now consider the expressions of the real and imaginary parts of A{ri.^ > ''es ) • The 
real part x{r^^ , r^^) satisfies 

^^ Z^-, (iV-l)7r\ 

a;(re2,re3) = -a+ (7 + d) I 1 - cos — 1 ^ -a + a^ 

hence, x does not depend on £2 nor on £3, and a; > for a < a*. 
On the other hand, the imaginary part y{rg^,ri,^) is 

y(re2 ,rej = (£27 + S3O) sm — 

thus y does not depend on a, and since sin - — j^ ''" > then, 

\y{r+,r-)\ <y{r+,r+) and y(r+,r+)>0. 

Finally, when a decreases from a*, then x increases from zero, and hence 4'{^) 
decreases from +00. The first value of y to satisfy (J16p will be y(r_|_,r+) since 
it has the largest absolute value. Hence the first pair of eigenvalues to cross the 
imaginary axis will be X(r^,r+)ei ~ '^(r_,r_)eii as required, while the real parts of 
all other eigenvalues, including A(r_^,r_)ei ~ ^(r^.r+)ei^ ^^^ ^till negative. D 

Note that from Lemmas [5] and |6l it follows that the first bifurcating eigenvalue 
for c > is precisely the non-resonant eigenvalue for c = 0, that has the largest 
imaginary part. 
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